Cooling a mechanical resonator by quantum interference in a triple quantum dot 
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We propose an approach to cool a mechanical resonator (MR) via quantum interference in a triple 
quantum dot (TQD) capacitively coupled to the MR. The TQD connected to three electrodes is 
an electronic analog of a three-level atom in A configuration. The electrons can tunnel from the 
left electrode into one of the two dots with lower-energy states, but can only tunnel out from the 
higher-energy state at the third dot to the right electrode. When the two lower-energy states are 
tuned to be degenerate, an electron in the TQD can be trapped in a superposition of the degenerate 
states called the dark state. This effect is caused by the destructive quantum interference between 
tunneling from the two lower-energy states to the higher-energy state. Under this condition, an 
electron in the dark state readily absorbs an energy quantum from the MR. Repeating this process, 
the MR can be cooled to its ground state. Moreover, we propose a scheme for verifying the cooling 
result by measuring the current spectrum of a charge detector adjacent to a double quantum dot 
coupled to the MR. 

PACS numbers: 03.65.Ta, 85.35. Be, 42.50.Gy 



I. INTRODUCTION 



Mechanical resonators (MRs) with a high resonant fre- 
quency and a small mass have wide applications and are 
attracting considerable recent attentions 0, [Ij. Tech- 
nically, these MRs can be used as ultrasensitive sen- 
sors in high- precision displacement measurements Q , de- 
tection of gravitational waves 3 or mass detection [jj. 
Also, quantized MRs can be useful in quantum infor- 
mation processing. Indeed, quantized motion of buck- 
ling nanoscale bars has been proposed for qubit implc- 
mentationlgl, @] and also for creating quantum entangle- 
ment |8Hl0j[. However, for all these applications, a basic 
prerequirement is that the dynamics of the MRs must 
approach the quantum regime. 

Quantum behaviors of a MR are usually suppressed by 
the coupling to its environment. One way to approach 
the quantum regime is to increase its resonant frequency 
so that an energy quantum of the MR becomes larger 
than the thermal energy. Recently, MRs based on metal- 
lic beams 11 1 and carbon nanotubes fl2| have been de- 
veloped, which have resonance frequencies of several hun- 
dred megahertz. However, for a MR with a frequency of 
200 MHz, a temperature lower than 10 mK (below the 
present dilution refrigerator temperature) is required to 
maintain the MR at the quantum regime. To attain the 
quantum regime, one needs to cool the MR further via 
coupling to an optical or an on-chip electronic system. 
Numerous experiments on cooling a single MR via radia- 
tion pressure or dynamical backaction have been reported 
(see, e.g., fl3l - l2"i| ). Theoretically, cooling by couplin g to 
a Cooper pair box or to a three-level flux qubit [23| 
via periodic resonant coupling have also been proposed. 
In these schemes, a strong resonant coupling between the 
MR and the qubit is required to cool the MR to its ground 
state. 



A. Sideband cooling of a MR 

In the weak coupling regime, a conventional method 
for cooling the MR is the sideband cooling approach (see, 
e.g., Ref. j24l43l| ). In this case, a MR is coupled to a two- 
level system (TLS) in which the two states can be elec- 
tronic states in quantum dots [2^ - |26| , photonic states 
in a cavit y |27H 30||. or charge states in superconduct- 
ing qubits }32j . In order to achieve ground-state cooling 
in the sideband cooling approach, the resolved-sideband 
cooling condition u) m 3> T (with u) m denoting the os- 
cillating frequency of the MR and F the decay rate of 
the TLS) must be followed in order to selectively drive 
the lowest sideband of the TLS. Then, the excitation of 
the TLS from the ground state \g) to the excited state 
|e) and the subsequent decay from this excited state to 
the ground state will, on average, decrease an energy 
quantum in the MR. This process can be described by 
\g, n) — > |e, n — 1) — > \g,n — 1), where n denotes to the 
state with n phonons. However, the frequency of a typ- 
ical MR is about 100 MHz [Til, O- It is in general of 
the same order of the decay rate of the two-level system. 
This indicates that the resolved-sideband cooling condi- 
tion is not easy to fulfill. Violating the condition means 
that the processes of a carrier transition and a subsequent 
sideband transition (i.e., \g,n) — > |e, n) — > \g,n + 1)) will 
occur. This will heat up the MR instead and suppress 
any ground-state cooling (29[. 



B. Cooling atomic motion via quantum 
interference in a three-level atom 

For laser cooling of atoms, an alternative approach [33| 
based on quantum interference in the internal degrees 
of freedom of the atoms without the need to follow the 
resolved-sideband cooling condition has been proposed. 



In this approach, an additional state is coupled to the ex- 
cited state 1 3) of the TLS to form a A-shaped three- level 
system [sec Fig.QJb)] . The two lower-energy states in this 
three-level system are tuned to be degenerate. Due to dis- 
sipation of the excited state, the atom will eventually ar- 
rive at a particular superposition of the two lower-energy 
states which is orthogonal to the excited state. This phe- 
nomenon results from the destructive quantum interfer- 
ence between the two transitions |1) — > |3) and |2) — » |3) 
and the superposition state is called the dark state j34| . 
When atomic motion is also considered, the carrier tran- 
sition (|1, n) — >• |3, n)) and thus the heating process of the 
atomic motion in the sideband cooling scheme is hence 
suppressed. It was shown that atomic motion can be 
cooled to its ground state in the non-resolved sideband 
regime [33j . 




C. Cooling a MR via quantum interference in a 
triple quantum dot 

In the present work, we propose a new scheme to cool 
a MR via capacitivc coupling to a triple quantum dot 
(TQD) schematically displayed in Fig.QJa). We consider 
the strong Coulomb-blockade regime so that at most one 
electron is allowed to present at one time in the TQD. 
The TQD acts as a three-level system in A configura- 
tion, in which the two dot states |1) and |2) (i.e., the 
single-electron orbital states in dots 1 and 2) are cou- 
pled to a third (excited) state |3) via two tunnel barri- 
ers [see Fig. Ufa)]. Here, the degrees of freedom of the 
MR is analogous to the motional degrees of freedom of 
the atoms discussed above and the TQD is an electronic 
analog of a three-level atom. We will show that by prop- 
erly tuning the gate voltages, one can degenerate the two 
lower-energy states and obtain a dark state in the TQD. 
By capacitively coupling to the TQD, the MR can be 
cooled to its ground state, in full analogy to the slowing 
down of the atoms via quantum interference. Compar- 
ing with the cooling of atoms [33{ , our approach has the 
following potential advantages: (i) Our cooling system is 
completely electronic and can be conveniently fabricated 
on a chip, (ii) Simply by adjusting the gate voltages, it 
is easy to achieve the two degenerate lower-energy states 
required for realizing destructive quantum interference in 
the TQD. Moreover, in contrast to the cooling of a MR 
by coupling it to a superconducting qubit [35| , the decay 
rate T of the higher-energy state of the TQD, which is 
equal to the rate of electrons tunneling from the TQD to 
the electrode, is tunable in our case by varying the gate 
voltage. 

Moreover, we also propose a method to verify whether 
the MR is successfully cooled by coupling it to a double 
quantum dot (DQD). This DQD can be reduced from 
the TQD by applying appropriate gate voltages. When 
the DQD and the MR is tuned into a strongly disper- 
sive regime in which the transition frequency difference 
between the two subsystems is much larger than the 



FIG. 1: (color online) Schematic diagram of a TQD system. 
Dots 1 and 2 are both tunnel-coupled to dot 3 (with interdot 
coupling strengths Sli and Q.2, respectively) while they are 
only capacitively coupled to each other. A MR is capacitively 
coupled to dots 1 and 3 of the TQD. (b) A three-level system 
in A configuration. The energy detunings between the two 
lower-energy states and the third excited state are respectively 
Ai and A2. The coupling strength between the state |1) (|2)) 
and the state |3) is Q,\ (Sl 2 ). 



coupling strength between them, the coupling between 
the MR and the DQD only yields a phonon-number- 
dependent Stark shift to the transition frequency of the 
DQD. This Stark shift corresponds to the shift of the 
resonant peak in the current spectrum of a charge detec- 
tor. Thus, by measuring the shift of the resonant peak, 
one can readout the phonon-number state and examine 
whether the MR is successfully cooled or not. 

This paper is organized as follows. Section II intro- 
duces a microscopic model for the coupled MR-TQD 
system. We show that the TQD is an electronic ana- 
log of a three-level atom driven by two electromagnetic 
fields. Also, we show how the TQD evolves into the dark 
state. In Sec. Ill, we derive a master equation to describe 
the quantum dynamics of the coupled MR-TQD system. 
With this master equation, we further derive in Sec. IV 
the master equation of the MR by eliminating the TQD 
degrees of freedom. Moreover, we calculate the steady- 
state average phonon occupancy of the MR and show 
that the MR can indeed be cooled to its ground state by 
using the quantum interference in the TQD. In Sec. V, 
we propose a method to verify if the MR is successfully 
cooled by measuring the full-frequency current spectrum 
of a charge detector. Section VI summarizes our results. 
In the appendix, we give a detailed derivation of the mas- 
ter equation for the reduced density matrix of the MR. 
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II. A MECHANICAL RESONATOR COUPLED 
TO A TRIPLE QUANTUM DOT 

A. Model 

The device layout of a MR coupled to a TQD is shown 
in Fig. [H[a). The TQD is connected to three electrodes 
via tunneling barriers. In the TQD, dots 1 and 2 are only 
capacitively coupled to each other and electrons cannot 
tunnel directly between them. Such capacitively coupled 
dots have already been achieved in experiments (see, e.g., 
[36|). In contrast, electrons can tunnel between dots 1 
and 3 as well as between dots 2 and 3. Here we focus 
on the strong Coulomb-blockade regime, so that at most 
a single electron is allowed in the TQD. Thus, only four 
electronic states need to be considered in the TQD, i.e., 
the vacuum state |0), and states |1), |2) and |3) corre- 
sponding to a single electron in the respective dot. The 
MR is capacitively coupled to dots 1 and 3 and this is 
schematically shown in Fig. QJa). 

The total Hamiltonian of the whole system reads 

-fftotal = H + Hi n t + Ht + H op . (1) 
The unperturbed Hamiltonian Hq is defined by 



Hq — -fflcads + Htqd + Hr + H p h, 



(2) 



where terms on the R.H.S. of Eq. (|2|) denote Hamiltoni- 
ans of the electrodes, the TQD, the MR and the thermal 
bath given by 



(3) 



ik 



Htqd = —A 1 a\a 1 -A 2 ala2 + (Qia\a 3 + fi 2 <4 a 3 + H.c), 

(4) 

H R = u m tfb, (5) 
H P h = ^Z u i b l b i- ( 6 ) 



We have put % = 1 and the energy of the state |3) is 
chosen as the zero-energy point. c\ k (cjfc) is the creation 
(annihilation) operator of an electron with momentum 
k in the ith electrode (i = 1,2 or 3) and a\ creates an 
electron in the ith dot. The phonon operators and 
b respectively create and annihilate an excitation of fre- 
quency uj m in the MR. In Eq. ([6]), the thermal bath is 
modeled as a bosonic bath with b q (b q ) being the creation 
(annihilation) operator at freqency uj q . 

The electromechanical coupling between the MR and 
dots 1 and 3 of the TQD is given by 



Hmt = -g (4 a 3 - a\ai){tf + b), 



(7) 



with a coupling strength g = rjuj m . For a ty pica l elec- 
tromechanical coupling, rj ~ 0.1 (see, e.g., Ref. j37|). The 
tunneling coupling between the TQD and the electrodes 



H T = ^2(n ik atctt+Kc.), 



(8) 



ik 



where £lik characterizes the coupling strength between 
the ith dot and the associated electrode via tunneling 
barrier. Moreover, the coupling of the MR to the outside 
thermal bath is characterized by 



H ep = J2n q (blb + R.c), 



(9) 



with a coupling strength tt q . 



B. Analogy between TQD and A-type three-level 
atom in two driving electromagnetical fields 

We now show that in the absence of the MR, our TQD 
system is analogous to a typical A-type three-level atom 
in the presence of two classical electromagnetical hclds. 
This hcld-drivcn three-level system is often used in quan- 
tum optics for producing a dark state (see, e.g., 34j|). 
The Hamiltonian of the field-driven A-type three-level 
system can be written as 

Ha = wiajai + L02a\a 2 + w 3 aja 3 
+fl a cos(us a t)(a\a,3 + ajai) 

cos(ojbt)(ala 3 + <4a 2 ), (10) 

where u>i (i = 1,2 or 3) is the energy of the ith state 
in the three-level system. Also, u) a and a;& are the fre- 
quencies of the two driving hclds and fl a and Of, are the 
corresponding driving strengths. In order to eliminate 
the time-dependence of the Hamiltonian in Eq. fj L0[) , we 
transform the system into a rotating frame defined by 



U R 



_ „iRt 



with 



R = oj a a\ai + Lo^a^a^ — u>3(a\ai + a\a,2 + ala^). (11) 



The transformed Hamiltonian is 



H A = U R 1 H A U R + iUn 1 U Rl 



(12) 



where the first term is evaluated as (within the rotating- 
wavc approximation) 

— l t t t 

U R Haur = LUia[ai + bJ2a 2 a,2 + ^30.30,3 

1 ^ Q 1 t , t i , ^ b i t 1 t \ 

+ -^-( a l a 3 + a 3 a l) + ^( a 2 a 3 + a 3 a 2j: 

(13) 

and the second term gives 

iUr X Ur = uj a a\ai + ujba\a,2 — u>3{a\ai + a\a 2 + a' 3 a 3 ). 

(14) 



Thus, Eq. ([T2J) reduces to 

H A = —Aia\ai - A 2 aja 2 

+Qi(a\a3 + a|ai) + Sl 2 (a2,a 3 + a\d2), (15) 

where Ai — u> 3 — u>i — to a and A 2 = u> 3 — w 2 — lui, are the 
frequency detunings while Oi = £7 a /2 and il 2 = 
are the effective driving strengths of the two fields. It is 
now clear that the Hamiltonian in Eq. (|T5|) is formally 
equivalent to that of the TQD in Eq. (QJ. This shows 
that the TQD we propose here is an electronic analog of 
a A-type three-level atom driven by two electromagnetic 
fields. 



C. "Dark state" in the TQD 

From the study of quantum optics, the existence of a 
dark state in a A-type three-level atom when the lower- 
energy states become degenerate, i.e., Ai = A 2 , is able to 
suppress absorption or emission. Below we demonstrate 
that a similar dark state also exists in the TQD [38| . 

After tracing over the degrees of freedom of the elec- 
trodes, the quantum dynamics of the TQD in the absence 
of the MR is described by 

pd = £tqdP<2 

= -i[H T QD, Pd] + riX>[aJ]p + T 2 V[a\]p d + T 3 V[a 3 }p dl 

(16) 

where pd is the reduced density matrix of the TQD and 
(i = 1, 2 or 3) is the rate for electrons tunneling into 
or out of the ith dot. The notation T> for any operator 
A is given by 

V[A]p = ApA*-^[A*Ap + pAlA]. (17) 

Considering equal energy detunings of the lower energy 
states |1) and |2) with respect to the excited state |3), 
i.e., Ai = A 2 = A, the eigenstates \g), |— ) and |+) of 
the TQD become 

\g) = p\3)-YiW) + Sh\m, 

h) = ^(fi 2 |l>-fii|2», 

|+) =a|3) + ^(fi 1 |l)+fi 2 |2)), (18) 

where a = cos(0/2), /3 = sin(0/2), tan 6 = 20/ A, and 
O = ^/f2 2 + f2 2 . The corresponding eigenenergies are 



E Z = 7T 

2 



= -A, E+ = -- 



(19) 



with </> — \J A 2 + 4£7 2 . For simplicity, we consider equal 
couplings of the three quantum dots to the correspond- 
ing electrodes, i.e., Y\ = T 2 = T 3 = V. Based on the 



Mi 




FIG. 2: (color online) Effective tunneling processes of elec- 
trons through a TQD represented in the eigenstate basis. 
An electron can tunnel from the left electrode into the three 
eigenstates |+), | — } and \g), with rates F/3 2 , F and Fa 2 , re- 
spectively. Note that the total tunneling rate is 2F because 
an electron tunnels from the left electrode to the TQD via 
two tunnel barriers (each having a tunneling rate T). In the 
eigenstate |+) (\g)), it will tunnel out to the right electrode 
with a rate To 2 (T/3 2 ). However, if the electron occupies the 
dark state |— ), no further tunneling occurs and the electron 
is trapped. 



eigenstate basis of the TQD in Eq. (TIB]) , the equations 
of motion for the reduced density matrix elements of the 
TQD are obtained from Eq. (fT6]l as 



poo = -2rpoo + Tfi 2 Pgg + Fa 2 p ++ + Ta(3(p +g + p g+ ), 
Pgg = Ta 2 poo - F/3 2 p gg - -a(3(p +g + p g+ ), 
P— = I> 



00: 



P++ = T/3 2 Poo - ra 2 p ++ - -aP(p +g + p g+ ), 

r 

P+g = ~i(E+ - E e )p +g - - P+g 

r 

- 2 Q ^++ + Pgg) - Ta PPoo- 



(20) 



Figure [2] schematically show effective electron tunneling 
processes through the TQD as described by Eq. (|20[) . 
Starting from an initially empty TQD, an electron can 
tunnel from the left electrode into any of the three eigen- 
states, with tunneling rates T/3 2 , T and Ta 2 for eigen- 
states |+), |— ) and \g), respectively. An electron in the 
eigenstate |+) (\g)) then tunnels out of the TQD to the 
right electrode with a rate Ta 2 (T/3 2 ). However, if the 
electron occupies the state |— ), no further tunneling oc- 
curs because, being orthogonal to |3), it is decoupled from 
the right electrode. Therefore, an electron in the TQD 
will be trapped in the state |— ), which is called the dark 
state in quantum optics (38| . This dark state results from 
the destructive quantum interference between the tran- 
sition |1) — > |3) (i.e., the electron tunneling from state 
|1) to state 1 3) in the TQD system) and the transition 
1 2) — i> 1 3) (i.e., the electron tunneling from state |2) to 
state 1 3)). 
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III. EFFECTIVE HAMILTONIAN AND 
MASTER EQUATION FOR THE COUPLED 
MR-TQD SYSTEM 

We now study the coupled MR-TQD system. Rather 
than analyzing directly the energy exchange between the 
MR and the TQD which involves tedious algebra, we 
apply a canonical transform U = e s on the whole system, 
where 



S = rj{a\a 3 — a\ai)(b — 
The transformed Hamiltonian is given by 



(21) 



= -ffioads + H ph + H cp + u m tfb - Aiajai - A 2 a|o 2 
+ [tt 1 a\a 3 B 2 + Q, 2 a\a 3 B + H.c], 
+ 2J [fii* a\ c lk B + D,2k 4 c 2fe 



-fl 3k a\ c 3 feB t +H.c], 



(22) 



where we have neglected a small level shift of ij 2 uj m to 
the states |1) and |3) and we have also defined 



B = exp[-r?(6 - & f )] 



(23) 



To describe the quantum dynamics of the coupled MR- 
TQD system, we derive a master equation (under the 
Born-Markov approximation) by tracing over the degrees 
of freedom of both the electrodes and the thermal bath. 
Up to second order in rj, the master equation can be 
written as 



dp 



-iu m [tfb, p] -i[H>TQT>,p] 

~i{V(tf ~b),p}+£ TP + £n P , (24) 



where 



V = f][2n 1 (a\a 3 - %ai) + ^ 2 (a 2 a 3 _ a l a 2)] , 



CtP 



£dP = 



TV[a\]p + TV[a\]p + TV[a 3 ]p 
+ ij 2 T 1 (V[a\tf]p + V[a\b]p) 
+ V 2 r 3 (V[a 3 tf}p + V[a 3 b}p), 
7[n(cj m ) + l]V[b]p + jn{u m )V[tf}p. 



(25) 



(26) 
(27) 



Here, the Liouvillian operator CtP accounts for the dissi- 
pation due to the electrodes and Cnp represents the dis- 
sipation at the MR induced by the thermal bath. Also, 7 
denotes the decay rate of excitations in the MR induced 
by the thermal bath and n(tu m ) is the average boson num- 
ber at frequency ui m in the thermal bath. 



IV. GROUND-STATE COOLING OF THE MR 



A. Master equation for the reduced density matrix 
of the MR 



In the limit 7 <C g <C uj m , the TQD is weakly coupled 
to the MR and can be regarded as part of the environ- 
ment experienced by the MR. The degrees of freedom of 
the TQD can then be adiabatically eliminated [13, 
and the master equation for the reduced density matrix 
p of the MR is given by (see Appendix) 

p, = -i{uo m + 8 m )[tfb, p} + 7^ {7H w m) + 1] + A-(u m )} 
•x[2bptf - (tfbp + ptfb)} 

+ \[ in {Lo m ) + A + {u m )][2tfpb - (btfp + pbtf)}, (28) 



where 5 m is the driving-induced shift of the MR fre- 
quency. In Eq. (|28p , the additional terms A + and A_ arc 
induced by the coupling with the TQD. With this mas- 
ter equation, one obtains the equation of motion for the 
phonon- number-probability distribution, p n = (n\p\n), 
of the MR: 



dp n 

dt 



{j[n(uj m ) + 1] + A_}[(n + l)p n+ i - np n ] 
+ [jn(uj rn ) + A + ][np n -i - (n + l)p„], (29) 



Moreover, the equation of motion for the average phonon 
number, (n) = ^2 n np n , in the MR can be obtained from 
Eq. ([29]) as 



dt 



■ = -(7 + W){n) +-/n(io m ) + A+, 



(30) 



where W = A_ — A + . In order to cool the MR, one needs 
W > (i.e., A_ > A+). 



B. Steady-state solution 

From Eq. (|30[) , the steady-state average phonon num- 
ber in the MR is 



"si 



771 (uj m ) + A 4 
7 + W 



(31) 



where the term 7n(w m ) in the numerator is due to the 
thermal bath while A + results from the scattering pro- 
cesses by the TQD. We assume that the MR is initially 
at equilibrium with the thermal bath, so that the initial 
phonon number in the MR is n(uj m ). In order to cool 
down the MR significantly, one needs a large cooling rate 
W ^> 7 to overcompensate for the heating effect of the 
thermal bath. At the end of Sec. IV, we show that this 
can be achieved using typical experimental parameters. 

Here we consider Ai = A 2 = A so that the dark state 
exists. The transition rates A± are found to be (see Ap- 
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pcndix) 



f> 2 4[o 2 - Wm K l ±A)] 2 + w , 2 „r 2 



(32) 

where p^ is the steady-state probability of an empty 
TQD. To cool the MR, one needs A- > A+, which is 
fulfilled cither when A > and < us m , or when A < 
and n > cj m . Assuming also W 3> 7, the steady- 
state average phonon number in the MR is approximately 
given by 



n st 



in{u> m ) 
W 



n f . 



Here n/ = A + /W which gives 

4[f2 2 - w m (ui m - A)] 2 



16Aw m « - ft 2 ) 

C. Optimal cooling condition 

It is easy to see that n / reaches the minimum 

r 

mm _ / J_\2 



(33) 



(34) 



(35) 



when the term in square brackets in the r.h.s. of Eq. (|34[) 
becomes zero, i.e., 



or 



SI 2 = u> m (uj rn - A), 



5 (A + . 



(36) 



(37) 



Therefore, by properly choosing the parameters f2, u) m , 
and A so that the optimal cooling condition in Eq. (|36|) 
is fulfilled and A 3> T, the steady-state average phonon 
number in the MR can be much smaller than unity, im- 
plying that ground-state cooling of the MR is possible. 
Moreover, the phonon number nj achievable according 
to Eqs. (p4|) and (|35|) is identical to the previous results 
for the cooling of trapped atoms via quantum interfer- 
ence [33|. However, the additional advantages of a solid- 
state cooling system proposed here are that it can be 
fabricated on a chip and is highly controllable. Specif- 
ically, all the relevant parameters (i.e., the detuning A, 
the tunneling rate V and the intcrdot coupling strengths 
ill and Q2) can be controlled by tuning the gate volt- 
ages in the TQD. Thus, for a fixed frequency uj m of the 
MR, the optimal cooling condition in Eq. (ptrJf can be 
conveniently fulfilled. 

The underlying physics of the optimal cooling condi- 
tion in Eq. (|36|) can be understood based on the eigen- 
state basis of the TQD. In the limit 7 <C g <C w m con- 
sidered here, the TQD arrives quickly at the dark state 
|-). The coupling between the MR and the TQD will 
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-4 -3 -2 

Energy detuning A/co m 

FIG. 3: (color online) Contour plot of the steady-state aver- 
age phonon number n Bt in the MR as a function of the nor- 
malized driving detuning A/o; m and the normalized interdot 
coupling 0,/ujm. The two solid curves correspond to n st = 0.05 
and 0.02. The black dasded line represents Q 2 = cu m (uj m — A) , 
under which the MR can be optimally cooled. We have chosen 
Q.i—Q.2= fi/v2 and typical parameters uj m = 27rxl00 MHz, 
r = u m , Q = 10'\ 77 = 0.1, and n(u) m ) = 21. 



excite the TQD to the state |+) most readily when the 
frequency uj m of the MR is equal to the transition fre- 
quency (<fi + A)/2 between the states |— ) and |+), i.e., 
w m = (cf> + A)/2. This corresponds to the transition 
\—,n) —> |+,n— 1). The excited electron subsequently 
tunnels to the right electrode, i.e., |+,n — 1) — > |0,n — 1). 
This whole process extracts an energy quantum from the 
MR. When this cycle repeats, i.e., |0, n) — > | — ,n) — > 
\+,n — 1) — > \0,n — 1) -> • • • , the MR is cooled to the 
ground state. Here we emphasize that the resonance con- 
dition for exciting the TQD from the state |— ) to the 
state I +) via the MR is equivalent to the optimal cooling 
condition in Eq. (|36|) . An electron can also relax from the 
dark state |— ) to the ground state \g) by releasing energy 
to the MR. However, this heating process of the MR is 
strongly suppressed because the frequency of the MR is 
off-resonant to the transition |— ) — > \g) in the TQD. 

Figure |3] displays a contour plot of the steady-state av- 
erage phonon number of the MR (n st ) as a function of 
the effective interdot coupling D. (= \JVt\ + S! 2 ) and the 
energy detuning A. Here we choose A < and > ui m 
to make sure that W > 0. For these typical parameters, 
a small n st < 0.05 is predicted over a wide range of values 
on the f2 — A plane. This implies that ground-state cool- 
ing of the MR should be experimentally accessible. Fur- 
thermore, to estimate the coo ling rate W, we use typical 
experimental parameters [IJ l40j : ui m — 27r x 100 MHz, 
A = -2ttx 300 MHz, and g = 2irx 10 MHz. The interdot 
couplings are chosen as Oi = — 27rx 141 MHz to fulfill 



the optimal cooling condition fi 



— A). Using 

Eq. (|32|) . one arrives at a cooling rate W w 2-7T x 2 MHz. 
Considering a MR with a quality factor Q = 10 5 (see, 
e.g., Ref. one has 7 = uj m /Q ~ 2n x 1kHz. There- 
fore, appreciable cooling with W 7 can be achieved. 



In this MR can be cooled from, e.g., an initial 

temperature T = 100 mK corresponding to n{io m ) =21 
down to T = 0.8 mK with n st = 0.017. 

In contrast, for sideband cooling of a MR [24T - I27T ] . the 
resolved-sideband cooling condition tu m 3> T must be fol- 
lowed for ground-state cooling of a MR. For a relatively 
large decay rate, only MR with a very high frequency 
(which becomes fragile in experiments) can be cooled. 
On the other hand, for cooling via quantum interfer- 
ence in the TQD proposed here, the cooling conditions 
O 2 = ui m (oj m — A) and A 2> T do not require a high MR 
frequency. 



V. A SCHEME FOR VERIFYING THE 
COOLING OF THE MR 

A. Quantum dynamics of coupled MR-DQD 
system in the presence of a charge detector 

To verify whether the MR is cooled or not, we pro- 
pose a scheme in which the MR is coupled to a two-level 
system realized by a double quantum dot (DQD). The 
state of the DQD is in turn measured by a nearby charge 
detector in the form of, e.g., a quantum point contact 
(QPC). This setup is schematically shown in Fig. 2) Ex- 
perimentally, the cooling and the verification setups are 
all fabricated on the same chip with shared components. 
Dots 1 and 3 in the TQD from the above cooling setup 
can make up the DQD after a gate voltage is applied to 
decouple dot 2 from dot 3. Also, the QPC should be 
decoupled from the rest of the system during cooling by 
applying a high gate voltage. 

The Hamiltonian of the system is given by 

H = H R + tf D QD + #QPC + #int + fldet- (38) 

The Hamiltonian Hr of the MR and the coupling H[ nt be- 
tween the MR and the DQD are already given in Eqs. (O 
and ©. Here i?DQD, #qpc and H^et are respectively the 
Hamiltonians of the DQD, the QPC and the coupling be- 
tween them and are given by 

A 

HuQD = — — Vz + filfx) 

#QPC = y^USkCs k C S k + UDqc]j a C Dq , 
kq 

H dct = £)(T - X^)(4 k c Dq + H.c), (39) 

where a z — 03123 — a\ai and a x = a^ai + a\a^ are the 
Pauli matrices. Also, Cik (cL) is the annihilation (cre- 
ation) operator for an electron with momentum k in ci- 
ther the source (i = S) or the drain (i = D) of the QPC. 
T is the transition amplitude of an isolated QPC and \ 
is the variation of the transition amplitude caused by the 
DQD. For simplicity, we assume that the dctunning A of 
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DQD QPC 




MR 

FIG. 4: (color online) Schematic diagram of a MR capaci- 
tively coupled to a DQD which is under measurement by a 
nearby QPC. The energy detuning between the two dot states 
in the DQD is A and the interdot coupling strength between 
them is fii. 



the DQD is zero. The DQD has the eigenstates 



|5) = ^-(|1>-|3», |e) = ^-(|l) + |3)), (40) 

with \g) (|e)) being the ground (excited) state. Rewriting 
the Hamiltonian in Eq. (|38|) on the eigenstate basis of 
the DQD, we have 

H = uj m tfb + Qip z - gp x {tf + b) + H QPC 

+ J2( T -XP x )(cl k c Dq + K.c), (41) 

kq 

where p z = \e)(e\ - \g)(g\ and p x = \e){g\ + \g){e\ are the 
Pauli matrices. 

We consider the coupled MR-DQD system in the 
strong dispersive regime where the coupling strength is 
much smaller than the difference between the transition 
frequency 2f2i of the DQD and that of the MR, i.e., 
g <C 5 = 20i — u rn . This regime was previously consid- 
ered to study whether the vibration of a MR coupled to a 
superconducting circuit is classical or quantum mechani- 
cal [4l| . In this regime, the phonon in the MR is only vir- 
tually exchanged between the DQD and the MR. Thus, 
the coupling of the DQD to the MR docs not change the 
occupation probability of the electron in the DQD, but 
only results in phonon-number-dependent Stark shifts on 
energy levels of the DQD. Moreover, the Stark shifts can 
be detected by measuring the full-frequency current spec- 
trum of the QPC. 

Applying both a rotating-wavc approximation and a 
canonical transformation U' = e s with 

s' =r)(p-tf - P+ b), n = g/S, (42) 
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to the Hamiltonian H, one obtains [up to 0(rj 2 



where 



H w uj, 



, t , , 2f2 1 + 5 2 (2&tb + l)/J 
,.0'oH / 



+y s (\9)(9\ + \e)(e\)+H QPC 



(43) 



From Eq. (|43|) . after taking the trace over the degrees of 
freedom of the QPC, one obtains the following master 
equation for the reduced density matrix elements of the 
coupled MR-DQD system [H: 



Pgn,gn 
Pen, en 

Pgn,en 



1+Pgn,gn + (7- + ld)P en. en 

1+Pgn,gn ~ (7- + id)p en, en 1 
~'iSfiPgfi,en ~ ("Tl ~o~) Pgn,ei 



(44) 



where S n = 2Qi + g 2 (2n + l)/5 and 71 = 2irg s gdX 2 eVd 
with g s (g^) being the density of states for electrons in 
the source (drain) of the QPC and Vd the bias voltage 
across the QPC. Here 7± = 7i(l=pA n ), with A„ = S n /eVd, 
are the QPC-induced excitation and relaxation rates be- 
tween the ground state and the excited state of the DQD. 
Also, 7^ is the relaxation rate resulting from the coupling 
of the DQD to the thermal bath. Since the dissipation 
rate of the MR is much smaller than that of the DQD, 
dissipation of the MR is neglected. In Eq. (|44p . the re- 
duced density matrix element pi n ,in (i = g, e) gives the 
occupation probability of the state \i,n) of the coupled 
MR-DQD system ,while Pi n j n (i ^ j) describes the co- 
herence between the states \i,n) and \j,n). The equa- 
tions of motion for other elements, e.g., pi n ,jn' ( n 7^ ti'), 
which are decoupled from those considered here, are not 
shown. Using Eq. (|44[) and the normalization condition 



Pn 



Pgn,gn\^) 



one finds 



(7- + ld)Pn 



270 
(7- + ld)Pn 
270 



Pgr: 



-270* 



Pen, en (0 ~. Pn \ ~ Pn Pen, en (0)1 ^ ^ 7 
270 27o 

Pgn,en\t ) — Pgn.en 



(45) 



where 70 = 7i + 7d /2 and p n is the probability that the 
MR is at state In). 



B. Current spectrum of the charge detector 

The dc current through the QPC is related to the elec- 
tron occupation probability in the DQD and is given 
by @ 



I{t) = eD Pll + eD'p 33 = \{D + £>') + |(£)' - D)(a z 



(46) 



D = 2ng s gd{T - xfVd, 
D' = 2^g s gd{T + xfV d , 



(47) 



are the respective rates of electron tunneling through the 
QPC when dot 3 is respectively occupied or empty [43| . 
Therefore, one can define the current operator as 



7(t) =I + ha z {t) =I + hp x (t), 



(48) 



with Jo,! = e(D ± D')/2. According to the Wiener- 
Khintchine theorem, the power spectrum of the current 
through the QPC is 0] 



S{u) = RcJ e^ T dT[(I(t)I{t + T )) - (7(t + r) (/(*))]. (49) 


Substituting Eqs. (jUJ) and (gSJ into Eq. (j49j), we get 



S(u)/S = 1 + 



27172 
7i + 72 
27i72 



^p„(l - np n ) 
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7i + 72 



K-Pn, 



To 



7 2 + (S n + L0f 



(50) 



where 72 = 2Trg s g d T 2 Vd, k = (7+ - t_ - j d )/2 r y , 
and Sq = 2eIo is the current-noise background. From 
Eq. (|5"0"|) . one sees that the current spectrum of the QPC 
consists of peaks at resonance points uj — ±5 n . These 
peaks have width 70 and heights increasing with the prob- 
ability p n . For instance, the peaks at the resonance point 



3 2 (2n+l) 



(51) 



is shifted by g 2 (2n + l)/5 from 20i. Thus, from this 
peak shift in the current spectrum, one can readout the 
phonon-number state of the MR. 

Figure 4 plots the current spectrum of the QPC with 
two different coupling strengths between the MR and 
the DQD. Results for three cases in which the MR is 
respectively in its ground state (n = 0.01 <C 1), the first- 
excited state (n = 1.0) or thcrmalized with an average 
phonon number n st = 1.0 are plotted. Each resonance 
peak in the current spectrum corresponds to a phonon- 
number state of the MR. For the thermally distributed 
case, the current spectrum shows several peaks, where 
the relative area under each peak gives the probability 
of the corresponding phonon-number state. As shown 
in Fig. SJa), the distance between two adjacent peaks is 
smaller than the intrinsic peak width in the weak dis- 
persive regime, i.e., 2g 2 /S < 70, and hence the measured 
spectrum shows an ensemble. In this case, the phonon- 
number state of the MR cannot be measured. In the 
strong limit (2g 2 / 5 > 70), however, the ensemble can be 
individually resolved [Fig. BJb)], which allows us to de- 
tect the phonon number and also to verify the cooling 
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FIG. 5: (color online) Power spectrum of the current through 
the QPC when the phonon number in the MR are respec- 
tively n — 0.01 (black solid line), n — 1 (red dashed line), 
or given by the thermal distribution (blue dotted line), i.e., 
p n = n" t /(l + n Bt ) n+1 with n st = 1. The coupling strength 
between the MR and the DQD is g = 0.1 cj m (a) and 
g — 0.3 u) m (b). The other parameters are ui m = 100 MHz, 
fii = 2cj m , 72 = 0.01cj m , 71 = O.272, and j d = 272. 



result of the MR. Indeed, a relatively strong coupling be- 
tween a MR and a quantum dot has been recently demon- 
strated [IH . The strong dispersive regime is thus achiev- 
able and one can apply the proposed coupled MR-DQD 
system to verify the cooling of the MR via measuring 
the current spectrum of a nearby charge detector (e.g., 
QPC). 



VI. DISCUSSION AND CONCLUSION 

Our proposal on ground state cooling of the MR re- 
quires that the TQD is able to evolve into the dark state. 
However, the dephasing of the TQD due to coupling to 
other degrees of freedom in the environment can project 
the TQD into one of the three localized states |1), |2), and 
1 3) and drive the system away from the dark state [38j . 
However, the dephasing between the two localized states 
|1) and 1 2) depends on the coupling strength and the en- 
ergy detunning between dots 1 and 2 [44j |. Here in our 
system no direct coupling exists between the two local- 
ized states |1) and |2), and thus the dephasing almost has 
negligible effects on the cooling efficiency of the MR. 

In summary, we have studied the cooling of a MR by 
capacitive coupling to a TQD. Wc show that when the 
two lower-energy localized states become degenerate, the 
TQD will be trapped in a dark state which is decou- 
pled from the excited state in the absence of the MR. 
With the MR in resonance with the transition between 
the dark state and the excited eigenstate in the TQD, 
we have shown that the MR can be cooled to its ground 
state in the non-resolved sideband cooling regime. More- 
over, we have proposed a coupled MR-DQD system in the 
strong dispersive regime for verifying the cooling result 



of the MR. In this regime, the coupling between the MR 
and the DQD induces a MR-phonon-number dependent 
shift of the transition frequency of the DQD. Thus the 
phonon-number state which characterizes the cooling re- 
sult of the MR can be detected by measuring the shifts of 
the resonance peaks in the current spectrum of a nearby 
charge detector. 
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Appendix A: Master equation for the reduced 
density matrix of the MR 

In this appendix, we derive the master equation [Eq. 
for the reduced density matrix of the MR from the 
master equation [Eq. (f2~4"|)] of the coupled MR-TQD sys- 
tem by eliminating the degrees of freedom of the TQD. 
In general, the dissipation rate of the MR is much smaller 
than the decay rate of the TQD, i.e., [n{uj m ) + l]j ^ T. 
The TQD hence attains its steady-state quickly and its 
perturbation to the MR can be regarded as part of the en- 
vironment [13, H3|. Up to the second order in n, Eq. (|24|) 
can be rewritten as 



t = Cp=[C ° 



Ci + £ 2 }p, 



(Al) 



where 



£oP = -iu) m [tfb, p] -i[H T q Bl p] 

+r 1 2?[a t 1 ]p + T 2 V[a\]p + T 3 V[a 3 }p, (A2) 
C lP = -i[V(tf-b),p}, (A3) 
C 2P = r?T x {V[a\tf\p + V[a\b}p) 

+„ 2 r 3 (V[a 3 tf}p + V[a 3 b]p) + C D p, (A4) 

are respectively the Liouvillians to zeroth, first, and sec- 
ond orders in ?y. At zeroth order in 77, the quantum dy- 
namics of the whole system is described by 



p{t) = £ oP (t). 



(A5) 



The MR and the TQD are decoupled. Since the TQD 
is at its steady state most of the time, one has p(i) = 
Pd ® Tr^{/o} with p s J" denoting the reduced density ma- 
trix of the TQD at steady-state and Tr^l- • • } the trace 
over the TQD's degrees of freedom. Equation (|A5|) has 
an infinite number of steady-state solutions. These so- 
lutions can be expanded in the basis of the eigenvectors 
n)(n'\ of the Liouville operator Cq with eigenval- 



Pd » I 
ucs X r 



-i(n — n')ojn 



CoPn 



An 



'Pn 
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Here, |n) (n = 0, 1, 2, ■ ■ ■ ) denotes the nth state of the 
MR and (n — n')ui m represents the energy difference be- 
tween the states |n) and |n'). For -q ^ 0, these states 
with different n are weakly coupled by the perturbative 
terms C\ and C 2 - To obtain the quantum dynamics of 
the MR, we project the system onto the subspace with a 
zero eigenvalue (n = n') of Cq. The projection operator 
V is defined by 



where 



C Vp = 0. 



(A6) 



Noting g <C u> m (i.e., i] <C 1), a second order perturbation 
expansion gives the following closed equation for Vp [25j ] 

00 

Vp(t) = VC 2 Vp(t) + J drVde^dVpit). (A7) 



Substituting Eqs. (|A3I) and (|A4|) into Eq. ([AT]) and taking 
the trace over the TQD degrees of freedom, the first term 
in Eq. (TAT]) becomes [H 

Tv d {V£ 2 Vp(t)} = l -[ in {u m ) +T?Ti$ \V[b\ii 

+\bWu m ) + 1] + r, 2 T iP f }V[b^, 

(A8) 

and the second term gives 



00 

Tr d | J dTVLie CoT LiP P {t)\ 



= —i(ui m + S m )[b'b,fx], 

+Re[G(iuj m )]V[b]p + Re[G(-i(j m )]V[b f ]fx, (A9) 

where p = TrjlVp} is the reduced density matrix of the 
MR and pff is the probability of an empty TQD at the 
steady state. Here, we have defined 

S m = lm[G(iuj m ) + G{-iw m )}. (A10) 

Thus, from Eqs. (fAT)) , (fA"H]) and (|A"9|) , one has 

P = -i{u) m + 8 m )[tfb,p] + -[jn(oj m ) + A + ]V[b]p 



+ -{l[n(u m ) + l]+A_}V[tf]p, 



A± = 2Re[G(±iw m )] + V^iPm- 



(A12) 



Eq. ()A11[) is simply the master equation (|28|) for the re- 
duced density matrix of the MR derived in Sec. IV. In 
Eq. (|A9|) . the correlation function G(s) is given by 



00 

G(s) = -Tr^y dtV(0)e c ^ nt V(0)( 



= -J dr(V(t)V(0))e st = -(V(s)V(0)), 


_ . (A13) 
where s = iuj m . Also, V(s) is the Laplace transform of 
V(t) and the Liouvillians £tqd is given in Eq. ()16|) . 

To determine the correlation function G(±s), one first 
calculate the Laplace transform V(s) of the interaction 
term V(t). For convenience, we introduce the vector op- 
erator a for the TQD whose components arc defined as 

a 1 = |l)(l|, 5- 2 = |2)<2[, a 3 -|3)(3|, 
*4 = |1><2|, 5-5 = |2><1|, 5-e = Il><3|, 
a 7 = |3)(l|, a 8 = |2)(3|, a 9 = |3}(2|, (A14) 

where the average value of each component is (<x;) = 
Tr{cTiPd}. Using this notation, one has 



V = 2r]Qi(& 7 - <j 6 ) + 77Sl 2 (og - <x 8 ) 
and thus 



(A15) 



G(s) = -2r ? fi 1 [5 7 (s) - S 6 (s)} - V n 2 [S 9 (s) - S s (s)}, 

(A16) 

where Si(s) = (ai(s)V(0)) with (ctj(s)) being the Laplace 
transform of (&i(t)}. From Eq. (TTo| . we find that (<7,(t)) 
obeys the equation of motion: 



dt 



= M(a(t))+B, 



(A17) 



(All) where 
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(A18) 



and B = (T 1 ,T 2 , 0, 0, 0, 0, 0, 0, 0) T . Here we have denned 
A d = A x - A 2 , Ax = -(iAi + |r 3 ), and A 2 = -(*A a + 
7^3). From Eq. (|A17[) . the steady-state solution of the 
vector (a) is calculated as 

(<7 st ) = M^B. (A19) 

Applying the Laplace transform to Eq. (|A17[) . one ob- 
tains 

s(a(s)) {*(&)) = M(a(s)) + -. (A20) 

s 

Moreover, the ith component of the vector (5(s)) is given 



by 

(a i (s)) = ^2L ik [{a k (0)) + ^}. (A21) 

Here the matrix L is defined as L = (si — M) _1 where I 
denotes the identity matrix. Assuming that the TQD has 
already attained its steady-state at initial time t = 0, i.e., 
(S-(O)) = (<7 st ), using Eqs. (1AT5|) . Sttfo . (IAT91) . (lASTI) 
and the quantum regression theorem [34| , one can obtain 
the correlation function G(s) and the scattering rates A± 
as given in Eq. (|A12[) . 
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